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Applying three independent techniques, we give numerical evidence for a finite superfiuid density in 
isotropic hole-doped t-J ladders: We show the existence of anomalous flux quantization, emphasising 
the contrasting behaviour to that found in the 'Luttinger liquid' regime stabilised at low electron 
densities; We consider the nature of the low-lying excitation modes, finding the 1-D analog of the 
superconducting state; And using a density matrix renormalization group approach, we find long 
range pairing correlations and exponentially decaying spin-spin correlations. 

PACS numbers: 74.72.-h, 71.27.+a, 71.55.-I 
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The behaviour of strongly correlated electrons con- 
fined to coupled chains is at present a topic undergo- 
ing much investigation; the reasons for this attention arc 
numerous. Firstly, with the behaviour of electrons un- 
der t-J or Hubbard type interactions in one-dimension 
now relatively well understood, the two-chain systems 
provide an interesting 'first step' towards the challenge 
of two-dimensions. Secondly, the unusual nature of the 
ground state of the undoped system, in particular the 
existence of a spin gaplil, leads to further interest with 
regards to 'gapped' superconducting behaviour. Further 
more, it is believed that compounds such as (V02)P207lI 
and SrCu203tJ may be described by a lattice of coupled 
chains. Whilst there is considerable literature concerning 
the possible phases in a t-J ladder, a complete picture 
is still far from being realised and our aim is to clarify 
the behaviour in a particular region, specifically consid- 
ering the nature of the gapped state when the system is 
doped. Various techniques have been applied previously 
in this hole doped regionB and there is some indication for 
hole pairing and modified d-wave superconducting corre- 
lations. 

The t-J Hamiltonian on the 2 x L ladder is defined as, 
H = J' J2(S jt i ■ Sj-,2 ~ jn jA n jt2 ) 
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where most notations are standard. (3 (=1,2) labels the 
two legs of the ladder (oriented along the a;-axis) while j 
is a rung index (j=l,...,L). We shall concentrate on the 
isotropic case where the intra-ladder (along x) couplings 
J and t are equal to the inter-ladder (along y) couplings 
J' and t'. 



At half filling, the hamiltonian reduces to the Heisen- 
berg modeLand the behaviour is generally relatively well 
understoodtJ. A simple interpretation is given by consid- 
ering the strong coupling limit (J = 0); in such a limit, 
the ground state consists of a singlet on each rung with 
a spin gap (~ J') which corresponds to forming a triplet 
on one of the rungs. With the introduction of interchain 
coupling J, the triplets can propagate and form a coher- 
ent band thereby reducing the spin gap. In the isotropic 
case the gap remains (~ 0.5J) and it is the nature of 
the state formed on doping such a system that we shall 
concentrate on. 

A possible phase diagram for the isotropic t-J ladder as 
a function of J/t and doping has been proposed recentlyD, 
the main features of which we summarize. Away from 
half filling, the spin gapped region persists, exhibiting 
hole pairing and, as we will show, possible superconduct- 
ing correlations. This behaviour is observed for a large 
region of J/t except for J/t >^ 2.5B where the system 
phase separates and also (perhaps) very small J/t where 
the gap may disappear. As the system is doped further, a 
Luttinger-like phase is stabilised exhibiting gapless spin 
and charge excitations. At very small electron densities, 
an electron paired phase exists. 

In this letter we will describe three independent forms 
of evidence for a finite superfiuid density in the spin 
gapped region of the phase diagram (we will work specif- 
ically with an electron density (n) = 0.8). The first set 
of results, and the ones we deal in most detail with, are 
based on the existence of anomalous flux quantization, a 
feature present in a superconducting state. Secondly, we 
consider the spin and charge excitation modes which may 
be used to characterize the state of a system. Finally, we 
present direct calculations of correlation functions which 
have been obtained using the density matrix renormal- 
ization group method. 

Our first set of results then concern the existence of 
anomalous flux quantization. The calculation involves 
threading the double chain ring with a flux $ and study- 
ing the functional form of the ground state energy with 
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respect to the threaded flux, namely Eq(&); throughout 
this letter we will measure the flux $ in units of the flux 
quantum <&o=/ic/e. In general, i?o ( < &) consists of an en- 
velope of a series of parabola, corresponding to the curves 
of individual many body states £L(<I>), exhibiting a pe- 
riodicity of one. Byers and YangQ have shown that in 
the thermodynamic limit, Eq(&) exhibits local minima 
at quantized values of <&, the separation of which is 1/n 
where n is the sum of charges in the basic group; these 
local minima in Eq($>) must be separated by a finite en- 
ergy barrier. Hence, for a superconductor we would ex- 
pect minima in Eq(<&) at intervals of 1/2; these minima 
are related to the existence of supercurrents which are 
trapped in the metastable states corresponding-, to the 
flux minima and are thus unable to decay awajcl. This 
phenomenon is known as anomalous flux quantization. 
It should be mentioned that the existence of anomalous 
flux quantization is an indication of pairing and is not 
sufficient in itself to imply a superconducting state. 

Detailed studies of the attractive Hubbard model on 
two-dimensional latticesa have indicated the presence of 
anomalous flux quantization, confirming the existence of 
superconducting correlations in the ground state; In con- 
trast, the repulsive Hubbard model exhibits no anoma- 
lous flux quantization. 

In addition to the existence of flux quantization, the 
function £'($) also gives a quantitative value of the su- 
perfluid density, defined in one dimension by 

A distinction should be made between LL (the super- 
fluid density) and D (the Drude weight )EHl: The super- 
fluid density corresponds to the curvature of the envelope 
of the individual many body states as a function of flux, 
whilst the Drude weight is obtained from the curvature 
of a single ground state many bod#-,energy level. In gen- 
eral these quantities are different^. In the thermody- 
namic limit no-particular applied flux is preferred when 
calculating DjEl and hence we consider the curvature of 
the whole LEq(<&) curve when considering the superfluid 
density. Note that the existence of superconductivity re- 
quires both anomalous flux quantization and a finite su- 
perfluid density (and infact D s has no real meaning in 
the absence of anomalous flux quantization). 

Numerically, the application of a flux through the dou- 
ble chain ring is achieved by modifying the kinetic term 
of the hamiltonian such that 

c ],P;s c 3+i,l3;s >-> c l/3;s c J+b/3;s ei2?_ (3) 

where $ is the flux through the ring measured in units of 
$o and L is the length of either chain. Hence the appli- 
cation of a flux is numerically equivalent to a change in 
the boundary conditions of the problem; $=0 represent- 
ing periodic boundary conditions and representing 
anti-periodic boundary conditions. 



The technique we have employed is exact diagonalisa- 
tion of finite systems, specifically 2x5 and 2 x 10 double 
chain rings with intermediate electron densities (n) = 0.8 
and m) = 0.4 corresponding to the regions of the phase 
diagram where we expect spin gapped or Luttinger liquid 
behaviour respectively. Note that the electron number is 
always a multiple of 4 in order to guarantee that antifer- 
romagnetic correlations are not frustrated when one goes 
around each chain. The modes of the system are char- 
acterized firstly by their spin: singlet and triplet excita- 
tions correspond to charge and spin modes respectively. 
It is also useful to consider the parity of the states of the 
system under a reflection in the symmetry axis of the 
ladder along the direction of the chains: Even (R x = 1) 
or odd (Rx — —1) excitations correspond to bonding 
or anti-bonding modes respectively. Finally, it is neces- 
sary to consider the momentum, kx=2im/L in order to 
determine the dispersion relation of each mode. Imple- 
mentation of these quantum numbers and symmetries is 
straightforward using exact diagonalisation methods and 
the various excitation modes may be obtained by calcu- 
lating the ground state energy in each symmetry sector. 

Concentrating initially on J/t— 0.5 (n) — 0.8 , we 
show in Figs. [j]a(b) all possible spin and charge modes of 
the 2 x 5 (2 x 10) system, for all possible momenta, as 
a function of applied flux. In the case of the larger sys- 
tem, we have omitted some of the details of the excited 
states to simplify the diagram, showing the full spectrum 
only for <!> < 0.25. For both system sizes, the minimum 
energy function Eo(&) is formed by charge (spin zero) 
bonding modes; the excited modes with different quan- 
tum numbers move further from the ground state as the 
system size increases (a result we have checked by finite 
size scaling techniques) and hence will not interfere with 
Eo(Q). The existence of minima at intervals of half a 
flux quantum (i.e. anomalous flux quantization) clearly 
indicates the existence of pairing. 

In order to probe the behaviour of Eq(&) further, we 
consider the quantity L [Eo($>) — Eq(& = 0)] as a func- 
tion of <J> for various values of J/t and (n) (L is the length 
of the ladder). Note that the curvature of this function 
in the thermodynamic limit gives the superfluid density. 
Fig. ^a shows the contrasting behaviour obtained when 
keeping J/t fixed at 1.0 and varying the electron filling, 
specifically (n) — 0.4 and (n) = 0.8 (both the 2x5 and 
2 x 10 results are shown). This plot clearly shows the 
existence of anomalous flux quantization for a filling of 
(n) = 0.8 and its absence for (n) = 0.4. The occurence of 
the absolute minima at different values of flux ($ = and 
$ = 1/2 for (n) = 0.8 and (rt) = 0.4 respectively) can 
be explained by considering the non-interacting Fermi 
sea for the two fillings; a lower energy state is formed 
by choosing the flux (and hence boundary conditions) to 
give a closed shell. Fig. ||d shows an equivalent plot but in 
this case keeping (n) constant at 0.8 and varying the pa- 
rameter J/t from 0.5 to 4.0. In this case anomalous flux 
quantization is exhibited for J/t— 0.5, whilst for J/t— 4.0 
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L \Eq{<&) — Eo(& = 0)] appears to scale to a flat function, 
consistent with the existence of a phase separated region 
(D s = 0). 

Except for the region believed to be phase separated, 
the form of the curve L [E (<&) — Eq ($ = 0)] appears to 
show only relatively small finite size effects and may be 
easily extrapolated to the thermodynamic limit; hence, 
by considering the curvature, an accurate value of the 
superfluid density (andJJie Drude weight) can be found. 
In a future publication^ we analyse the specific values 
in more detail but for this letter we emphasise that D s 
scales to a finite value in the thermodynamic limit in the 
regions which are not phase separated. 

The second form of evidence for a superfluid density 
lies in the dispersion behaviour of the low energy spin and 
charge excitations of the system. The possible phases of 
a particular model can be characterized by the number of 
charge and spin modes which are gapless at zero momen- 
tum. The one-dimensional analog of a superconductor 
has one gapless charge mode, a gap to all spin excitations 
and dominant pairing. The phase diagram for the various 
possible phases for a chain with a Hubbaxd (rather than 
t-J) hamiltonian has recently been foundtj and interest- 
ingly on doping away from half filling, this spin gapped 
phase with one gapless charge mode is stabilised, exactly 
the situation we will suggest for the t-J case. In addi-i 
tion, both this article and a recent paper by NagaosatLa 
note that in this gapped phase, in contrast to a Luttinger 
liquid phase, the power law term in the density-density 
correlation function at 2k f is missing. 

In Fig. U we show the dispersion of all the spin and 
charge modes (for both bonding and antibonding sym- 
metry sectors), for the 2 x 10 ladder with (n) = 0.8, 
J/t = 1.0, corresponding to the region of the phase dia- 
gram where we believe the superfluid to exist. There are 
several obvious features of this diagram. Firstly there 
is a finite gap to spin excitations; secondly there is (at 
least) one vanishing charge mode (bonding) as k x >-^0 and 
thirdly, there is no sign of 2k f charge gapless modes. 
All of these features point towards a spin liquid with 
dominant superconducting correlations as previously ex- 
plained. We should note the possibility of a new charge 
gapless mode at finite momenta corresponding to the 
fluctuations of the pair density as indicated by the dip 
of the bonding charge mode in Figure ||. Some more de- 
tailed descriptions of these results (along with the corre- 
sponding results from other regionSfOf the phase diagram) 
are given in a separate publication^. 

The final results we present are direct numerical calcu- 
lations of various correlation functions obtained using the 
density matrix renormalization group approachlla; this 
technique allows much larger systems to be considered 
than is possible with exact diagonalization techniques. 
We present three types of correlations: Firstly the equal- 
time rung-rung pair field correlation function ^A^A^ 



indices and 1 and 2 indicate the chain; this correlation 
function then creates a singlet pair on rung j and removes 
a singlet pair from rung i, a direct measure of the motion 
of hole pairs in the spin gapped state. The second corre- 
lation function that we consider is the spin-spin correla- 
tion function • S| along one of the chains. Finally we 
measure the density-density correlations along one of the 

chains defined by (pi,iPj,i) where p iyl = Y,a c i,i;a c i,U°- 

In Fig. [| we show the results of a direct calculation of 
these correlation functions for a 2 x 30 ladder with open 
boundary conditions and (n) = 0.8 and J/t = 1.0 i.e. 
corresponding to the region where we expect superfluid 
behaviour. The data is shown on a log-log plot and the 
dashed line corresonds to a slope of —2, the dotted to 
a slope of —1. The most immediate feature of the re- 
sults is the fact that the pairing correlations are much 
longer range than the other correlation functions, decay- 
ing with behaviour slower than 1/L. The density-density 
charge correlations appear to decay close to 1/L 2 whilst 
the spin-spin correlations are consistent with an expo- 
nential decay (as we would expect for the spin gapped 
state). These results are then further evidence for our 
interpretation of the behaviour in this region of the phase 
diagram: dominant superconducting correlations with a 
spin singlet (gapped) wavefunction. 

In summary, we have presented three independent 
forms of evidence for a superfluid density in hole doped t- 
J ladders. Firstly we have shown the existence of anoma- 
lous flux quantization and a well-converged and finite D s ; 
Secondly we have studied the low lying modes, finding a 
spin-gapped state with a gapless charge mode and no 
gapless 2k f excitations; Finally we have presented direct 
calculations of correlation functions, showing long range 
pairing correlations and exponentially decaying spin cor- 
relations. 
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FIG. 1. Energy as a function of flux (in units of &o=hc/e) 
for different system sizes with J/t= 0.5,(n) = 0.8. We show 
all possible momenta for various quantum numbers: For the 
charge modes the solid lines correspond to bonding and the 
dotted lines to anti-bonding, whilst for the spin modes the 
dashed lines correspond to bonding and the dot-dashed lines 
to anti-bonding. Figure la (lb) corresponds to a system size 
of 2 x 5 (2 x 10); for the larger system size we give only the 
charge bonding modes and the lowest lying spin anti-bonding 
mode in full in order to simplify the diagram. 

FIG. 2. L[E (<f>) - E (<f> = 0)] where L is the length of 
the ladder and Eo(<&) is the ground state energy with an ap- 
plied flux <E>. The dashed lines correspond to 2 x 5, the solid 
lines to 2 x 10. Fig.2a shows the results for (nj — 0.4 and 
(•n} = 0.8 both with J/t = 1.0, whilst fig. 2b shows the results 
for J/t = 0.5 and J/t = 4.0 both with (n) = 0.8. 

FIG. 3. Spin and charge excitation modes of a 2 x 10 ladder 
versus momentum k x (in units of 7r); (n} = 0.8 and J/t = 1.0. 
The quantum numbers associated with the various symbols 
are shown on the plots. 

FIG. 4. Log-log plot of various correlation functions versus 
| i — j | (real space separation) for a 2 x 30 open chain with 
(rij = 0.8 and J/t = 1.0. The dashed line has a slope —2 and 
the dotted line —1. The correlation function are explained in 
the main text. 
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a) 2x5 J/t=0.5 <n>=0.8 
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